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A  MATIIl'MATICAL  MODIX  OI«  Till*  PROBABILITY  OK  PHRFORATION 
Ol'  Tin:  HUMAN  SKULL  BY  A  BALLISTIC  PROJIlCTILK 


1,  INTRODUCTION 

In  civilian,  as  well  as  niiliUiry  environinenis,  there  are  soinetiines  hazards  I'rom  ballistic 
iVagments.  Because  of  the  c\lremc  vulnerability  of  IliC  brain  to  penetrating  injury,  the  dilTerencc 
between  a  relatively  harmless  wound  and  a  lethal  one  often  hinges  on  whether  the  skull  is  perforated 
wlicn  struck  by  one  of  these  projectiles.  In  the  military  application,  an  important  part  of  determining 
the  vulnerability  of  the  soldier  or,  conversely,  the  effectiveness  of  protective  helmets  lies  in  our  ability 
to  predict  the  circumstances  under  which  the  skull  will  withstand  the  impact  of  a  fragment  or  bullet, 
In  this  report  a  mathematical  model  is  developed  which  provides  an  esCmate  ot  the  probability  of 
perforation  of  the  inner  table  of  the  skull  and  subsequent  penetration  of  the  brain  by  a  cube  or 
sphere.  This  probability  is  considered  i  function  of  the  striking  velocity  and  other  physical  properties 
of  the  projectile  as  well  as  the  thickness  of  the  bone  under  the  impact  area. 

DI  RIVATION  OF  THH  MODEL 

The  skull  is  composed  of  inner  and  outer  surfaces  of  hard  compact  bone  and  an  area 
between,  called  cancellous  bone,  which  is  less  dense  because  it  contains  a  higher  percentage  of  living 
cells.  The  impacting  projectile,  when  it  perforates,  punches  out  a  hole  in  the  outer  table  roughly  its 
own  size.  The  area  over  which  the  force  is  exerted  spreads  within  the  cancellous  portion,  causing  a 
larger  and  more  irregular  hole  in  the  inner  table.  This  phenomenon  is  often  called  “cratering”  and  is 
used  by  forensic  pathologists  to  tell  exit  from  entrance  wounds  in  the  skull.  This  is  possible  because 
the  larger,  irregular  hole  is  in  the  outer  table  with  an  exit  wound. 

To  simplify  the  mathematical  model,  this  phenomenon  is  idealized  by  the  following 

assumptions: 


1 .  In  the  region  of  the  impact  the  skull  is  a  flat  homogeneous  plate  of  bone  with 
uniform  mechanical  properties. 

2.  When  the  skull  is  perforated  by  a  projectile,  a  truncated  right  circular  cone  of  bone 
is  driven  out  in  front  of  the  fragment. 

3.  The  top  of  the  truncated  cone  has  an  area  equal  to  the  mean  presented  area  of 
the  projectile:  i.e.,  the  projectile  orientation  is  assumed  to  be  random. 

4.  The  velocity  of  sound  in  bone  is  much  greater  than  the  velocity  of  perforation,  so 
compressional  effects  may  be  neglected. 

5.  The  probability  of  perforation  depends  only  on  the  logarithm  of  the  ratio  of  the 
stress  imposed  on  the  skull  by  the  impacting  projectile  to  the  tensile  strength  of  the  bone. 

6.  The  underlying  probability  distribution  governing  perforation  is  well  approximated  by 
the  logistic  distribution  on  the  model  variable,  the  log  stress  ratio  of  assumption  5  (e.g.,  if  the  true 
distribution  were  Gaussian  this  assumption  would  be  valid).  Some  of  the  factors  contributing  to  the 
variance  of  the  probability  distribution  are  variations  in  the  mechanical  properties  of  bone  from  skull 
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to  skull,  dilTorcuccs  in  the  proportion  of  cortical  to  cancellous  bone,  and  random  deviations  I'rom 
some  or  the  assumptions  listed  above  (for  instance,  deviation  I'rom  conical  shape  in  the  crater),  each 
of  which  is  unique. 


A  schematic  diagram  of  the  penetralin(>  projectile  is  shown  in  lipure  I .  l-or  purposes  of 
illustration,  the  projectile  is  idealized  as  a  rinht  circular  cylinder  impacting  eiul-r)n.  The  terms  which 
will  be  used  in  developinji  the  model,  most  of  which  appear  in  I'itiure  1,  are  dcllned  in  the  list  of 
symbols.  The  units  used  are  not  SI  but  aiv  those  customarily  useil  in  wound  ballistics  and  arc 
retained  here  as  a  matter  of  convenience. 


I 

I 


Figure  1.  Idealized  Schematic  Diagram  of 
Cylinder  Perforating  Skull 

The  surface  area  of  the  side  of  the  truncated  cone  is  given  by: 

S  =  rr  (r  +  R)s 

But  s  =  t  sec^  and  R  =  r  +  t  tan'P,  so 

3  =  ir  (2r  +  t  tan^)t  sec^' 


(1) 


C) 
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Tlic  minimum  averiipe  force  cxcrtcti  in  pushing;  out  the  conical  section  of  skull  may  be  approximated 
by  dividint:  the  kinetic  energy  of  the  projectile  by  the  hypothetical  elastic  limit  of  the  bone,  i.e.,  the 
distance  that  the  plug  could  be  moved  before  it  v^mthl  ba*ak  loose  from  the  surrounding  tnaterial. 
Thus, 

force  ’s  '/2inv“/d 


tensile  strc'ss  ^  force/S 


stress  ratio  =  tensile  stress/tensiie  strength  ( 1) 


_ '/j  mv^/d 

T  TT  (2r  +  t  tan)t  scc'l' 


(3) 


file  quantity  2r  may  be  renlaced  by  the  variable  I)  =  2y'A/7r  to  generalize  to  a  projectile  with  a 
nonciicular  presented  area.  When  inserted  into  the  probabilitv  function,  the  multiplicative  et)nstanls  of 
equatii'ti  3  nui\  be  absorbed  into  logistic  parameter  a  (see  ctiuation  5  belov.).  With  these  modifications 
tlie  stress  ratio  becomes  the  model  variable, 

'/i  mv“  '/:  mv“/t“  * 

^  t~sec^(D/t  +  tan'I')  sec'KD/t  +  tan'l') 

We  will  use  the  logistic  probability  distribution  function  to  estimate  the  probability  of  penetration  of 
the  skull. 


P  = 


I _ 

I  4-  g-(a  +  blnX) 


(5l 


where  a  and  b  are  determined  from  data  (see  next  section). 

The  natural  logarithm  (In)  of  X  is  used  instead  of  X  itself  to  equalize  the  variance  over  the  domain  of 
definition.  This  requirement  is  most  easily  explained  by  use  of  an  example:  Suppose  we  are  concerned 
with  The  probability  that  a  certian  projectile  will  perforate  some  material  in  sheets  of  different 
thicknesses  where  the  probability  of  perforation  is  a  function  of  velocity  alone.  On  a  thin  sheet  where 
the  V5Q,  the  velocity  at  which  the  probability  of  perforation  is  0.5,  is  lOO  m/s.  we  would  expect  a 
standard  deviation  on  the  order  of  tens  of  meters  per  second.  On  a  much  thicker  sheet,  where  the 
Vi;Q  is  1 000  m/s,  we  no  longer  expect  the  standard  deviation  to  be  measured  in  tens,  but  in  hundreds 
of  meters  per  second.  In  other  words,  we  expect  the  standard  deviation  to  be  roughly  proportional 
to  the  magnitude  of  the  mean.  In  this  case,  the  variance  may  be  approximately  equalized  over  the 
whole  domain  by  dealing  with  the  logarithms  of  the  numbers  rather  than  with  the  numbers  themselves. 


*  The  equation  is  expressed  in  this  manner  as  an  aid  to  plotting  X  at  a  later  stage. 


II 
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The  data  used  were  vrhlaiiied  Irom  (wo  previrnis  studies  on  skull  penetration.*,**  Since 
ilte  analysis  on  these  studies  showed  that  the  response  of  dried  skulls  was  sijinirieantly  dillerent  from 
tl'.at  of  Iresh  skullcaps,  only  the  latter  were  used  in  detenninintt  the  parameters  u  and  h  lor  the  lojjistie 
luiK'tion  of  equation  The  physical  properties  of  the  Hve  projectiles  used  are  listed  in  the  table. 


Table.  l*roiee(iles  Studied  (All  Maile  td  Steel) 


Project  ilo 

1 

H  -  ■ 

Mass 

Mean 

dimension* 

Mean 

pivsenleil  urea** 

gm 

cm 

0.8.x  er,iin  sphere 

0.05-S 

().2.^8 

0.0445 

2  i  grain  cube 

0.135 

0.265 

0  1050 

4.2  grain  cube 

0.274 

0.333 

0,1 060 

lo  grain  cube 

1 .029 

0.514 

0.3966 

225  erain  cube 

14.694 

1 .2.36 

2.2933 

*  Oiaiuetcr  for  sphere;  edge  for  cubes. 
**  Surface  arca/4;  see  appendix  A. 


Tlie  ability  of  the  model  to  properly  scale  the  probability  of  perforation  of  skulls  of 
different  tliiekness  by  projectiles  of  different  mass  was  tested  against  the  of  the  above  five  pro¬ 
jectiles.  Tlte  v^Q  was  obtained  by  r'-eraging  the  five  highest  nonperforating  velocities  and  the  five 
Knu'st  perforating  velocities  for  each  skull  thickness  for  which  there  were  sidficicnl  data.+  Tlrese  data 
arc  sliowM  plotted  in  terms  of  the  model  in  figure  2.  Note  that  the  abscissa  is  the  logarithm  of  the 
numerator  of  the  model  variable  of  equation  4,  while  the  ordinate  is  the  logarithm  of  the  denom¬ 
inator.  A  straight  line  of  slope  1  represents  a  constant  ratio  of  the  numerator  and  denominator  tsee 
appendix  B).  When  drawn  through  the  data  in  figure  2,  this  line  allows  a  visual  assessment  of  the 
eoodness  of  fit  of  the  data  from  low  energy  to  high  energy.  The  energies  arc  mv-^Q  for  all  v.^q 
points.  Since  'I'  is  a  variable,  it  may  be  adjusted  to  cause  the  sphere  and  cube  data  to  be  on  a 
common  line.  A  'I'  of  23°  for  the  sphere  and  30°  for  the  cubes  not  only  superimposed  the  two 
groups  hut  tilted  the  individual  data  sets  to  form  the  best  straight  line  of  slope  1 .  Since  the  co¬ 
linearity  of  the  data  with  the  reference  line  is  relatively  insensitive  to  'T,  the  difference  (V°  between 
spheres  and  cubes)  is  more  significant  than  the  magnitude  of  4'. 


•  Millci.  J.  F.,  Ashman,  W,  P,,  and  Jameson,  J.  W.  Edgewoud  Arsenal  Technical  Report  FAIR  4.t73.  Ballistic  Limits 
of  Skulls  Against  Steel  Cubes.  April  1970, 

**  Mickiewicz,  A.  P.  Chemical  Systems  Laboratory.  Private  communication. 

T  To  obtain  the  maximum  number  of  Vjq  points,  a  few  were  averages  of  4  and  4  instead  of  5  and 


SKULL  BALLISTIC  LIMIT 
X  =  .85  SPHERE  V  ^  23  deg 

*  =  2.  I  CUBE  l|l  =  30  deg 

+  =4.2  CUBE 


Art  .  the  iletcrminiitioa  of  all  the  variables  in  etination  4  are  known.  It  is  then 
posMblo  to  lit  the  ilichotomous  ilata  (perl'oration  or  nonperloratitjn)  to  the  logistic  function  by  the 
..pproximate  least  squares  methoil  of  Walker  anil  Duncan.*  In  this  techni(|ue  all  data  are  used  to  approx¬ 
imate  the  lor  each  skull  thickncss/projcctile  combination.  I’or  each  shot  the  inavs,  striking 
velocity.  4'.  skull  thickness,  and  mean  presented  area  are  used  to  calculate  X  Irom  equation  4.  A 
second  number  is  associated  with  this  value:  I  lor  perl'oration.  0  lor  nonix'rf'oratMin.  Tlie  Walkcr-Duncan 
method  then  iteratively  converges  on  values  of  a  and  b  which  appioximalcly  miniini/e  the  sum  of 
squared  vliU'erenees  between  the  predicted  probability  of  penetration  (equation  5)  and  the  assigned 
N.diie  tU  to  h  lor  all  shots  by  all  projectiles.  The  resulting  values  ol  a  and  b  are: 

a  -  -K^.{)35.  b  =  4.7').^ 

lines  ol  slope  I  in  I'igure  2  represent  constant  values  ol  X  and.  therefore,  discrete  levels 
of  probability  of  perforation.  To  find  the  line  of  probability  we  insert  I*  =  0.5  m  equation  5 

iiul  solve. 


I 

p  =  -  =  Q  5 

I  +  e~^“  ^ 


or 


^-(a  +  blnXjQ)  =  I 


rins  implies 


a  +  bInXjQ  =  0 
Vi  mv~5Q 

InXcQ  =  In  ; - — - 

t-sec'l'(D/t  +  tan^') 

where  X^q  is  from  equation  4  with  v  =  vjq. 

Tlierefore. 


-a/b 


(6) 


ln[t“secW(D/t  +  tan^*))  =  a/b  +  Injli  niv^Q)  ("') 

It  IS  this  line  which  is  drawn  in  figure  2.  The  good  agreement  between  the  line,  derived  from  all  data, 
and  the  plotted  points,  derived  from  a  few  selected  data  points,  lend  some  credence  to  the  validity  of 
using  those  approximate  vjq  points  in  determining  the  value  of  'I'.  An  example  of  the  use  of  the 
model  is  given  in  appendix  C. 


*  Walker,  S.  It.,  and  Duncan,  D.  B.  Estimation  of  the  Probability  of  an  Event  as  a  Function  of  Several  Independent 
Variables.  Biometrika  54.  167-179  (1967). 
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4. 


SUMMARY  AND  CONC  LUSIONS 


T!k'  moilol  ilciivod  ahovc  is  based  on  the  iissiinipti«)n  that  a  pint!  ol  honionencous  material 
111  the  sliapc  ol'  a  truncated  cone  is  driven  out  of  the  skull  by  the  perl'oratinn  projectile.  'I'lie  area  ol' 
tile  small  oiul  of  the  truncateil  cone  is  assumed  equal  to  the  mean  presented  area  ol  the  projectile. 
Under  these  assumptions,  both  spheres  and  cubes  may  be  accommodated  in  the  same  model  by  letting 
the  hall  anelc  of  the  cone  take  on  din'eren.  values  I'or  the  two  shapes,  lixperieiice  with  previous 
models  ol'  penetration  has  shown  that  cluinkv.  irrei'ular  I'ratiments  usually  behave  in  a  manner  more 
>imilai  to  cubes  than  spheres.  I  herel'ore,  it  .,  v  .j»};ested  that  the  cube  hall'  anple,  'l^  =  .10',  be  used 
'or  I'laements  as  weii. 

riie  ilata  used  in  I'ittiuj;  this  motlel  were  obtaineil  I'rom  impacts  on  bare  bone.  In  predicting 
the  probability  ol'  penetratini’  the  skull  ol'  live  humans,  this  moilel  neglects  the  protection  ofjered  by 
the  scalp  and  (sometimes)  hair.  Thus  the  model  is  conservative  in  the  .sense  ol  predicting  somewhat 
higher  probabilities  ot'  penetration  that  would  actually  be  observed  in  live  human  skulls  ol  the  same 
thickness, 


C  aution  must  be  exercised  when  applying  these  criteria  to  projectiles  larger  than  those  in  the 
data  base  or  ot'  density  much  dit't'erent  from  that  of  steel.  There  are  no  data  available  to  test  the 
model  outside  the  range  of  masses  and  presented  areas  of  the  projectiles  of  the  table.  Large  projectiles 
at  lower  velocities  are  particularly  risky  because  of  significant  curvature  of  the  skull  over  large  areas 
and  different  mechanisms  of  skull  fracture  which  occur  under  these  conditions.  It  should  also  be  noted 
that  assumption  3  requires  random  orientation  of  the  impacting  projectile.  This  will  be  valid  for 
projectiles  whose  mean  dimension  is  not  much  larger  than  the  thickness  of  the  skull,  since  in  per- 
lorating  the  skull  the  projectile  will  ultimately  present  its  mean  area,e.g.,by  rotation  during  penetration. 
For  projectiles  whose  mean  diniension  is  greater  than  the  skull  thickness  (for  example,  paper  weights), 
t!ic  projectile  is  unlikely  to  perforate  presenting  its  mean  area,  and  the  effect  of  orientation  of  the 
projectile  will  be  important. 

5.  RECOMMENDATIONS  FOR  FURTHER  WORK 


When  examining  the  skullcaps  perforated  by  cubes,  one  does  n^t  see  flat-sided  holes  with 
nci’t  square  corners,  but  rounded  holes.  This  indicates  that  the  assumption  of  equality  between  the 
mean  presented  area  of  the  projectile  and  the  area  of  the  hole  in  the  outer  table  should  instead  be 
an  assumption  of  proportionality.  Differences  between  spheres  and  cubes  could  then  be  accounted  for 
by  different  constants  of  proportionality  rather  than  different  half  angles  Irregular  fragments  would 
have  a  proportionality  constant  somewhere  between  those  of  the  cube  and  sphere.  Unfortunately,  there 
are  no  data  available  for  irregular  fragments  nor  even  sufficient  data  for  the  spheres  and  cubes  to 
allow  an  accurate  determination  of  the  proportionality  constants. 


If  additional  data  are  acquired,  particularly  fragment  data,  it  is  suggested  that  the  following 
model  be  tried: 


mv^ 

^  ^  t(kD  +  t  tanW) 


where  k  depends  on  the  sh^pe  of  the  projectile  (sphere,  cube,  chunky  fragment,  etc.)  and  ^  is  a 
constant,  independent  of  f'^ape. 


APIM-NDIX  A 


Mi-AN  i»ri;si:nti;i)  ari:a 

Lot  us  first  ilovolop  tho  "simple  case”  of  the  mean  presented  area  of  a  sphere.  We  want  to 
express  this  as  a  fraction  of  its  total  surface  area.  To  obtain  the  surface  area  we  will  integrate  in  polar 
co-ordinates  as  in  llgure  A-l.  Notice  that  the  width  of  the  element  of  area  dA  is  r  sin  0  d^.  The  sine 
Ainction  is  necessary  because,  like  the  outer  surface  of  an  orange  segment,  it  must  narrow  to  a  point 
at  the  top  (0  =  0).  I'ven  if  the  wedge  is  infinitesimally  thin,  it  still  must  be  shaped  like  a  wedge. 
Another  way  of  looking  at  this  effect  is  to  imagine  0  held  constant;  then  wlien  we  rotate  around  the 
vertical  (/.)  axis  through  an  angle  d^i,  the  distance  moved  on  the  surface  is  r  sin  0  d^.  The  surface  area 
of  tile  sphere  is 


sin0  d0  d^  =  r^ 


f-cosO)  1^  d^4 


Referring  again  to  figure  A-l,  we  now  wish  to  obtain  the  mean  area  projected  on  the 
plane  PI  at  the  right  of  the  figure,  parallel  to  the  x-z  plane.  The  width  and  height  of  dA  as  projected 
on  this  plane  are  r  sinfl  sin^  d«i  and  r  sin0  d0.  Thus  the  projected  area  is 


Ap 


sin\i  sin^0  d0  d^ 


(A-2) 


where  the  limits  of  integration  on  ^  are  0  and  ir  because  the  region  from  tr  to  2n  is  on  the  far  side 
of  the  sphere  from  the  plane  PI  and  does  not  contribute  to  the  presented  area  required.  Thus. 


The  mean  presented  area  Ap  (as  expected)  is  the  same  as  the  cross-sectional  area  and  is  a  quarter  of 
the  surface  area.  The  angles  6  and  (i  may  be  interchanged  without  affecting  the  results. 
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y 


Figure  A-1.  Element  of  Area  of  a  Sphere  and 
Its  Projection  on  a  Plane 


We  now  turn  to  the  problem  of  the  mean  presented  area  of  a  polyhedron  with  an  arbitrary 
number  q  of  plane  elements  as  its  surface.  Let  the  area  of  a  typical  element  be  E.  Then  the  projected 
area  of  tile  element  E,  under  all  possible  orientations  of  the  polyhedron  with  respect  to  the  plane,  of 
projection  (x-7L  is  equivalent  to  the  projected  area  of  E  rotated  about  its  center  of  mass  as  in 
figure  A-2.  The  co-ordinate  system  is  centered  at  the  center  of  mass  in  such  a  way  that  the  normal 
n  to  the  area  E  makes  an  angle  9  with  the  z-axis.  Tlie  projection  of  the  normal  onto  the  x-y  plane 
makes  an  angle  (fi  with  the  x-axis.  The  integral  of  the  area  projected  onto  PI  is 


rf\ 

Jo  Jo 


sin0  sini^  sin0  d0 


(A-4) 


TT 

sin^  sin^0  d0  d<i 


ttE 


where  the  first  sin0  is,  as  above,  necessary  to  provide  equal  weight  to  all  possible  orientations.  Because 
I-  is  a  finite  area  (not  an  infinitesimal),  the  area  Aj  is  an  integral  projected  area  -  not  a  mean.  To 
obtain  the  mean,  we  must  divide  Aj  by  the  total  solid  angle  subtended  by  all  possible  orientations 
of  the  element  E  (i.e.,  the  solid  an^e  swept  out  by  n),  namely  4?:.  Tlius  the  mean  projected  area 
Ap  =  A,/47r  =  E/4. 


Appendix  A 
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Figure  A-2.  Element  of  Area  of  a  Polyhedron 
and  Its  Projection  on  a  Plane 


If  there  are  q  elemental  areas  in  the  polyhedron  each  having  a  mean  projected  area  equal 
to  a  quarter  of  its  surface  area,  the  total  mean  projected  area  is  a  quarter  of  the  total  surface  area, 
i.e.. 


q 

=  L 

i=l 


Aj/4 


=  -  E 

4 


-  4  As 


(A-5) 


Tliis  formula  is  valid  for  convex  surfaces  only.  If  there  were  a  concave  portion  on  the 
surface,  it  would  be  shielded  by  another  part  of  the  surface  on  part  of  the  interval  O-tt  of  one  or 
both  angles  of  equation  A-4.  Any  convex  shape  can  be  approximated  to  any  specified  accuracy  by 
making  q  sufficiently  large  and  the  Aj  sufficiently  small.  Thus  the  mean  presented  area  of  a  convex 

1 

solid  is  —  of  its  total  surface  area. 

4 


0 
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APIM'NDIX  B 


COMMENTS  ON  Tllli  SLOIM-  OJ-  A  UX;  -  LOO  PIOT 

Consider  the  graph  ip  the  figure,  where  the  two  variables  N  and  I)  are  plotted  on 
logarithmic  axes.  Suppose  the  data  indicate  a  straight  line  relationship  on  this  graph. 

Then. 

In  I)  =  p  In  N  +  Cl 
=  In  NP  +  Cl 
.-.  In  NP  -  In  D  =  -<i 


..  after  exponentiating. 


D 


If  p  =  1,  there  is  a  constant  ratio  between  N  and  D, 

If  p  =  2,  there  is  a  constant  ratio  between  and  D. 

If  N  =  'i:mv-/t“  and  D  =  sec'I'(D/t  +  tan^)  the  N/D  =  X  (see  equation  4).  A  line  of  slope  1  in  the 
plot  of  In  N  V  In  D  then  represents  a  constant  ratio  between  Vimv^/t^  and  sec'KD/t  +  tan^r). 


Figure.  General  Linear  Relationship  Between  In  N  and  In  D 


AIMMtNDIX  (’ 


i;xAMi»u;  oi<  Till-;  usi-  or  Tin-  matiii-matk  al  modi-l 

Lot  US  assume  that  wo  wish  to  know  the  velocity  corrospouilinn  to  a  0.5  probability  ol 
ponotratiou  lor  a  0.5*jiram  cube  <>!'  mean  presenteil  area  0.245  cin“  impacting  a  skull  of  0.5-cni 
tliioknoss.  From  oiiuation  b  we  obtain 

=  |2t“  sec'l'(l)/t  +  tan'l')exp(-a/b)/m  I 


whore 


I)  =  2  y/Kiii  =  2  v'0.245/ir  =  0.5585  cm 


=  30° 


t  =  0.5  cm 

m  =  0.0005  kg 

a  =  -13.035 


b  =  4.793 


Substitution  above  gives  V5Q  =  172  m/s.  Suppose  we  wish  to  know  the  probability  of 
penetrating  the  skull  if  the  cube  of  our  example  impacts  at  200  m/s.  From  equation  4  we  obtain 


X  = 


Vi  mv^ _ 

t^sec^(D/t  +  tan^) 


20.44 


Tlien 


1 

p  - - - - -  0  807 

1  +  exp(-a  -  blnX) 

Tlie  following  table  can  be  constructed  using  similar  calculations,  assuming  the  same  fragment  and 
skull  thickness; 


Velocity  m/s 


Probability  of  penetration 


150 

175 

200 

225 


0.209 

0.537 

0.807 

0.973 


The  value  of  probability  of  penetration  above  does  not  indicate  the  proportion  of  the  skull  thickness 
penetrated  but  rather  the  proportion  of  hits  which  would  be  expected  to  punch  out  a  plug  of  bone 


trcMij  the  skull  given  a  large  number  of  hits.  Thus,  in  the  example  above,  at  v  =  200  m/s,  80  7%  of 
all  shots  fired  would  be  expected  to  penetrate  for  that  fragment /skull  thickness  combination,  while 
I  would  not.  No  attempt  is  made  to  .-redict  the  outcome  of  individual  shots 
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